function [elements,processInfo] = InitializeSolutionStep(elements,model,processInfo)
ne = model.countElements();
dim = 2;
for i=1:ne
    e = model.getElement(i-1);
    f = e.getFace();
    ue = f.getUHatLocal();
    nn = f.countModes();
    ng = e.getNumGaussPoint();
    B = zeros(3,dim*nn);
    for j=1:ng
        xi = e.getGaussPoint(j-1);
        B = e.computeB(f,B,xi);
        C = elements(i).constitutiveMatrix(:,:,j);
        %compute stress
        elements(i).stress(:,j) = C*B*ue;
    end
end
end
